Instability of coherent states of a real scalar field 
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We investigate stability of both localized time-periodic coherent states (pulsons) and uniformly 
distributed coherent states (oscillating condensate) of a real scalar field satisfying the Klein-Gordon 
equation with a logarithmic nonlinearity. The linear analysis of time-dependent parts of perturba- 
tions leads to the Hill equation with a singular coefficient. To evaluate the characteristic exponent 
we extend the Lindemann-Stieltjes method, usually applied to the Mathieu and Lame equations, 
to the case that the periodic coefficient in the general Hill equation is an unbounded function of 
time. As a result, we derive the formula for the characteristic exponent and calculate the stability- 
instability chart. Then we analyze the spatial structure of the perturbations. Using these results 
we show that the pulsons of any amplitudes, remaining well-localized objects, lose their coherence 
with time. This means that, strictly speaking, all pulsons of the model considered are unstable. 
Nevertheless, for the nodeless pulsons the rate of the coherence breaking in narrow ranges of ampli- 
tudes is found to be very small, so that such pulsons can be long-lived. Further, we use the obtaned 
stability-instability chart to examine the AfHeck-Dine type condensate. We conclude the oscillating 
condensate can decay into an ensemble of the nodeless pulsons. 

PACS numbers: 03.65.Pm, 05.45.Yv, ll.lO.Lm, 11.27.+d 



I. INTRODUCTION 

Nonlinear localized field configurations, solitons, are 
currently considered as models of various physical ob- 
jects, from elementary particles and collective excitations 
in condensed matter to giant lumps of dark matter in the 
form of soliton stars and galactic halos 0,110. Stability 
properties of solitons were investigated by many authors, 
and a number of important results has been obtained 
(see, e.g., and references therein). In particular, Ho- 
bart [3 and Derrick @ have proved that static multidi- 
mensional scalar solitons are energetically unstable, and, 
hence, these objects cannot last in a real world for a 
long time. One way to avoid this theorem is to invoke 
time dependence. Along this line main efforts were fo- 
cused on the stability analysis of the stationary states, 
i.e., coherent states of a complex scalar field oscillating 
harmonically in time. It turned out, however, that for a 
wide class of relativistic models these states can be only 
conditionally stable, i.e., stable with respect to a certain 
type of perturbations (e.g., conserving the scalar charge) 
01 . As to the time-periodic states of a more general form, 
both complex and real, there are presently no strong an- 
alytical results on their stability. 

In this paper we examine stability of time-periodic con- 
figurations of the form 

= 0o(i,r) =a(t)u(r) (1) 

satisfying the nonlinear Klein-Gordon equation 

A0 + C/'(0) = O. (2) 

These are coherent states in the sense the field oscillates 
synchronously at all spatial points. It is necessary to 
stress that we consider real solutions, so that the energy 
density oscillates as well (in contrast to the stationary 



states for which a{t) oc e*'^*). Solitons with oscillating 
energy density are usually called pulsons. 

It turns out that for real (j) the ansatz ^ determines 
uniquely the potential U{4>) in Eq. Namely, if nei- 
ther a{t) nor u{r) are constants, the only potential ad- 
mitting such solutions will have the form (3 



(3) 



where rri^ and A are arbitrary constants. 

Originally, the Klein-Gordon equation with a loga- 
rithmic potential of this type has been introduced in 
the quantum field theory by G. Rosen |^. Later on 
Bialynicki-Birula and Mycielski have rediscovered this 
equation and also considered its nonrelativistic version, 
the nonlinear Schrodinger equation [Tol |. 

In inflationary cosmology and in modern supersym- 
metric field theories the logarithmic nonlinearities ap- 
pear naturally when quantum corrections to effective po- 
tentials are allowed for 11^, 13. 14\ . In this context 
the expression in the square brackets of Eq. ||2Jl can be 
treated as a dynamic inflaton mass term nig that is the 
bare inflaton mass term plus the logarithmic correc- 
tion. It can be represented in the commonly considered 
form [11 111 111 IH 111 E m EH by the substitution 
ln(j)'^ = 1 + ln($/M)^, A/m^ = —K, where $ is an infla- 
ton scalar field, M is a large mass scale, K is a. constant 
(usually negative and small). Thus our consideration is 
also relevant to dynamics of the pulson exitations of a 
real infiaton field oscillating around a vacuum value. 

Note that the multidimensional pulsons probably exist 
in other scalar models as well. Thus the long-living oscil- 
lating spherically symmetric localized states were numer- 
ical ly f ound in the sine-Gordon, (j)*, and (j)^ — (f)"^ models 
El El El El El (see 01 for a review). Unfortunately, 
the analytic form of these solutions is so far unknown. 



2 



The model (0) and Q is unique in the sense it has a 
whole family of exact pulson solutions of the form , all 
existing in any number of spatial dimensions H. This is 
also true for complex version of the model 29] . The 
real pulsons we are dealing with arc the limiting states of 
the complex ones, when the scalar charge tends to zero. 
Other limiting states are Q-balls for which a{t) oc e"^* 
[Tsl . It is believed that Q-balls can arise due to frag- 
mentation of the Affleck-Dine condensate [T l ll9ll20ll2l f . 
We will see below that the parametric instability of the 
oscillating condensate leads to the resonant fragmenta- 
tion that can give rise to the pulson formation at the 
nonlinear stage. Like Q-balls pulsons interact elas- 
tically or inelastically in collisions depending on their rel- 
ative velocities, phases, and rest masses 30, ol]. Thus, in 
model (I2Jl and Q the light pulsons with given relative ve- 
locities interact always elastically, independently of their 
phases. In contrast, the collisions of heavy pulsons can 
result in formation of the so-called explosons, localized 
states with exponentially growing amplitude [33|. For 
the intermediate masses the picture depends essentially 
on the phases of the colliding pulsons and im pac t velocity 
determining the duration of the interaction |31j . 

The above results suggest that there is a domain of 
parameters where pulsons are stable, at least in short 
time interactions. But in what sense? How long a pul- 
son conserves its characteristic features once interaction 
ends? If pulsons are long-lived objects they will be inter- 
esting candidates for the dark matter constituents having 
time-dependent density. What is known about stability 
of an isolated pulson at the long time scale? Surpris- 
ingly, but very few. In Ref. it was argued in favour 
of its perfect stability. No deviations from the exact solu- 
tion were found after about one thousand oscillations. 
However, our preliminary numerical experiments (32t i33l] 
have shown that the pulsons of certain amplitudes, even 
perturbed by computer round-off errors only, gradually 
lose their coherency, remaining well-localized oscillating 
objects. This has motivated the closer examination. 

In the present paper we clarify how long the pulsons 
can conserve the coherency depending on their param- 
eters. For this purpose we investigate stability of the 
spherically symmetric pulson solutions ^ with respect 
to small initial perturbations of an arbitrary form. 

The paper is organized as follows. In Sec. II the main 
properties of the real pulsons of the model considered 
are reviewed. Section III is wholly devoted to the linear 
stability analysis. We arrive at the singular Hill equa- 
tion and generalize the Lindemann-Stieltjes method to 
evaluate the characteristic exponent. On this basis we 
examine stability of the pulsons and discuss fragmenta- 
tion of the oscillating AfFleck-Dine type condensate. In 
Sec. IV we make some remarks concerning the complex 
pulsons and summarize the main results. 



II. PULSONS AS COHERENT STATES 

Assuming A positive, let us first eliminate the constants 
and A from consideration by the scaling t — + X^^^^t, 

In the new variables the 



A 



-1/2, 



)exp 



2A 



field (j) may be thought of as satisfying Eq. |(2Jl with the 
potential 



(4) 



It is the potential we will deal with. It has local minimum 
at (/) = and two maxima at (/i = ±1, at the minimum 
the potential having the singularity: its second derivative 
tends to infinity as ~> 0. 

The substitution of the ansatz into Eq. (jSJ leads 
then to two independent equations. 



d 
da 

d 

du 



a2(l -lna2) 



-u2(lnu2 _ 1) 



(5) 



(6) 



Note that the potentials in the square brackets of Eqs. 
(0) and jnj have the same form as the potential taken 
with plus and minus signs, respectively. The existence 
of the oscillating localized solutions Q is thus apparent 
from consideration of motion of a mechanical particle in 
these potentials. 

Let us consider in more detail the oscillatory solutions 
of Eq. (|SJ). Using the Hamiltonian and denoting ^ = 
(0 < flinax < 1, — 1 C ^ 1)) we obtain 

*2 



e'lnf. 



where 



(7) 



(8) 



In the case of small amplitudes, a: 
gives 

^(t) « COStiJo*- 



2 

max 



< 1, > 1, Eq. 



(9) 



Thus, we have quasi-harmonic high-frequency oscillations 
which are however nonlinear since their period, 



2n 



llnOmax 



(10) 



depends on the amplitude [32 
tion from Eq. ((TJ) we find 



|l/2' 

In the next approxima- 



T : 



2tt 



0.307 
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1 



(11) 



1, 



In the case of near-critical amplitudes, when a^^. 
the oscillations become almost rectangular and have the 
period 



T w 2\/21n- 



1 



(12) 
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FIG. 1: Oscillatory solutions of Eq. ^ for flmax = 0.1, 

Smax = 0.7, and Omax = 0.9999. 




FIG. 2: Spectrum of the spherically symmetric TV-nodal so- 
lutions «iv(r) of Eq. ©: uo(0) = e^/^ ui(0) = 9.726, 
U2{0) = 15.084, U3(0) = 20.526, .... 



Examples of solutions of Eq. (0) are shown in Fig. 1. 

The spatial structure of a pulson is determined by Eq. 
©. In the spherically symmetric case this equation has 
a discrete spectrum of localized A'^-nodal solutions MAr(r) 
with the first derivatives vanishing at the origin 

M (see 

Fig. 2). The simplest of them, the nodeless solution, has 
a Gaussian-like shape, 

«o(r) = e(3-'-^)/2, (13) 

and is usually called gausson |^, 0| . It is agreed that 
its effective radius equals v^- In the multinodal solu- 
tions, as r increases, the field undergoes spatial oscilla- 
tions of the half-wavelength L < 2\/2 and then decays 
as 

-(r-p„)V2 



{r:^r~pr^, iV > 1), (14) 



where Cat is the value of the last extremum of UAr(r) 
attained at r = p^, Cn (— l)^e^/^ {N oo), pj^ ~ 
NL. Thus, the pulsons of the model l(2Jl and 10} are 
well-localized states of the inhomogeneity length L, at 
all points the field oscillating coherently with the period 
T. (To return to the physical units these scales should be 
multiplied by A^^^^.) In our dimensionless variables the 
pulsons are characterized by two parameters only: the 
amplitude amax and the number of the nodes N (or T 
and UAr(O), respectively). 

It should be stressed that, due to nonanalyticity of 
U {4>) dX 4> — 0, the right-hand sides of Eqs. © and © 
are nonanalytic when a and u become zero. Hence, the 
solutions a{t) and u[r) themselves become nonanalytic 
at those points t and r where they pass through zero. 
Thus in the solution ^ we have dropped the terms which 
are small (of the order of ijOq"^) but nonanalytic when 
^(t) = 0. In general case from Eq. 10 it follows that ^(t) 
passes through zero at i = t,n as 



(15) 



where ± sign is taken for Ct(im) ^ 0. It is seen that ^^jj 
becomes infinite as t — > Similarly, one can show that 
in the vicinity of the n-th node {n — 1 , . . . , N) 



u(r) 



Ur{rn){r - Tn 

1 



1 



6 

+o(C 



[r - r„) ln(r - r„)^ 



(16) 



As we will see below, the nonanalyticity of U{(t)) gives 
rise to some specific features of the stability analysis. 



III. THE LINEAR STABILITY ANALYSIS 

Consider a small fluctuation iri{t, r) around the spheri- 
cally symmetric pulson JQ, = 0Q(t, r) -I- ri{t, r). In the 
linear approximation the equation for r/ reads 



(17) 



Seeking a solution in the form ri{t,r) cc X{t)^>{r) we 
arrive at the equations 



Xtt + {E-2-\iia'^)X = 0, 
A* + {E + lnu"^)^ = 0, 



(18) 
(19) 



where E is some constant. 

The expression in the brackets of Eq. H17|l is —U"{4)q). 
It becomes infinite, as well as the expressions in the 
brackets of Eqs. ifTFI) and lfTO|) . at the points tm and 
r„ where a{t) and u{r) become zero. Thus we need to 
analyze the second order differential equations with sin- 
gular coefficients. We begin with Eq. (|18|l which has the 
periodic singular coefficient In (t) and hence belongs to 
the class of Hill equations. 
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A. Singular Hill equation and generalized 
Lindemann-Stieltjes method 

It turns out to be very useful to look at the problem as 
a whole, considering first the Hill equation of a general 
form 

Xu + h{zit))X ^0. (20) 

We will assume that h{z) is an integral function of z, 
while z{t) is a real- valued periodic (of a period r) even 
fimction of t, having, in general, singularities, but such, 
that h{z{t)) remains still integrable. 

It is well known that the Hill equation describes the 
physical systems in which the parametric resonance can 
occur. In the context of our stability analysis we will 
be interested in real resonant solutions of Eg. 112011 . In 
accordance with the Floquet theory (see, e.g., [S^l), any 
one of these solutions can be represented as a linear com- 
bination of the fundamental solutions 

X+{t) = ^{t)e^\ X^{t)^^{~t)e-^\ (21) 

where Lpit) is a r-periodic or r-antiperiodic real function, 
/i > is the characteristic exponent. In the case that 
z{t) is unbounded it is impossible to obtain the solutions 
and evaluate /i by expansions in Fourier series, following 
the standard Hill approach. Another way is to apply the 
Lindemann-Stieltjes method j36i]. In some cases it allows 
one to obtain the results in a closed analytical form [s^, 
IsM I39I l40 l| . We first outline this method in the context 
of the general Hill equation H20|l with an extension to 
the case that the periodic function z{t) is unbounded. 
In doing so we follow the paper where the method 
was used to construct the resonant solutions of the Lame 
equation. 

The main idea is as follows. Let us treate 2: as a new 
"time" variable instead of In each interval of mono- 
tonicity of z{t) we define 

y{z)^X{t). (22) 

Assume that the periodic function z{t) satisfies the 
equation 

z^t^9{z), (23) 

where g{z) is an integral function of z. Eq. H2()|l then 
becomes 

9{z)y" + \g'{z)y' + h{z)y = (24) 

(hereinafter the prime denotes d/dz). 

Let us first suppose z[t) is bounded. Equation 
then shows that it is differentiable. Zeros of the function 
g{z) on the complex z plane, taken to be isolated, are 
singular points of Eq. (f^Hl . Since z{t) is periodic and 
real-valued, among singular points there are two, Ci and 
C,2i lying on the real axis and being minimal and maximal 



values which z{t) acquires at the end points of the inter- 
vals of monotonicity. Also, it follows that g'(Ci 2) 7^ 0- 
Physically, this is well understood, since C^i and ^2 can 
be treated as turning points in periodic motion of a me- 
chanical particle, e.g., of a nonlinear oscillator, under the 
action of the force g' /2. From Eq. (|23|) it is clear that the 
interval [Qi, ^2] does not contain other singular points of 
Eq. (123 . 

For example, in the case of the Mathieu equation we 
have z{t) = cos^ i, g{z) — 4z(l — z), so that Eq. has 
the regular singular points z ^ Ci =0, z — C,2 — both 
being the turning points. In addition, the equation has an 
irregular singularity at infinity. For the Lame equation 
z(t) = ST?{t,>c), g{z) = 4z(l — z)(l - )<^z). Equation 
H24II then has the regular singular points z = Ci = 0, 
z — C,2 — z = >c~'^ > 1, first two of them being the 
turning points, and a regular singularity at infinity. 

In general, it is easy to verify that the turning points 

2 are regular singular points of Eq. (|24|l , the exponents 
at each being and 1/2. This implies that in the vicin- 
ity of each turning point C, there exist two independent 
solutions of Eq. y^°\z-X) and y(^/^^(z;C), having 

asymptotics I + 0{z - Q and (z - + 0{z - ()], 

correspondingly. 

Now let us consider any one interval of monotonicity 
of the T-periodic even function z{t). Denote as yi{z) and 
2/2(2) those two linearly independent solutions of Eq. H24|) 
one of which coinsides, by Eq. H22I) . with the increasing 
solution (|21|l . X+{t)^ and another with the decreasing 
one, X^{t), on the interval chosen. Since ip{t) is either 
T-periodic or r-antiperiodic, the product ip{t)ip(—t) = 
X^X- — 2/12/2 = w(z) is always r-periodic even function 
defined on the whole t axis. Hence, at the end points 
of the intervals of monotonicity of z(t), i.e., at tm = 
mT/2 (m = 0,±1, . . . ), the derivative [{X^X^)t]t=trr^ = 
or, what is the same, 

i^'V9).=C^.. = 0. (25) 

In the vicinity of a turning point C the solutions 2/1 and 2/2 
can be represented as linear combinations of the solutions 
2/^°^ and y'^^^'^\ Consequently, the singularity [z — C)^^^ 
is the only one which the function w(z) = 2/12/2 might 
have. But its existence is in contradiction with Eq. 125|1 . 
because g'[C) ^ and, hence, g{z)z^c^ ^ 9'{C){z — C)- 
Therefore, the product 2/12/2 is analytic at z = 2- Re- 
call now that the interval [C1JC2] does not contain other 
singular points of Eq. and the singular points are 

assumed to be isolated. We thus conclude that on the 
complex z plane there exists a vicinity of the interval 
[C1JC2]! an open domain D D [C17C2]; which 2/12/2 
is an analytic function of z. In addition, it follows that yf 
and j/2 of necessity have singularities of the type (z— C)^/^ 
and, thus, cannot satisfy Eq. (|25|1 . 

Now consider the case that one of the turning points 
or the both are at infinity. This implies that at the corre- 
sponding instants tm the functions z{t) and Zt{t) become 
unbounded, the latter changing the sign. Nevertheless, 
we assume that in the vicinities of tm the function h{z(t)) 
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in Eq. (|20|l is integrable and X is continuous, whence 
it follows that Xt and, therefore, Wt are also continu- 
ous. Hence, as before, {■wt)t=t,„ = due to evenness and 
periodicity, so that we arrive at Eq. 125|) again, where 
Ci = — oo and/or C2 = +00. 

It is easy to verify that the bilinear combinations yf, 
yiy2, and ?/| constitute the fundamental system of solu- 
tions of the third-order differential equation 

g{z)w"'+y{z)w"+l^-g"{z)^Ah{z)^ w'+2h'{z)w - 0. 

Equation (|25|l is thus a common criterion for selection of 
the solution 



w = yiy2 



(27) 



from the set of solutions of Eq. (|26|l . In the case that z{t) 
is bounded, Eq. H25|l is the equivalent to the requirement 
that a solution of Eq. H26|) be analytic in D. If z{t) is 
unbounded, Eq. will give the boundary conditions 
at infinity which must be satisfied in solving Eq. H26|l . 
In this case w{z) will be analytic in a vicinity D of one 
of the intervals (—00,^2]) [CijOo), (—00,00). 

Thus, in a neighbourhood of any one point C G Z) we 
can write the expansions 




0^ 



(28) 



Substitution of H28|l into Eq. H2t)|l leads to the following 
set of equations for the coefficients: 



m+2 



m 



^ n{m -f n)g„„„+2W„ + 4 ^ (m + n) Wn = 

(29) 



n— 1 n— 

(to = 1,2,...). Thus for to = 1 we have 



Qgows + 3giW2 + (52 + 4:ho)wi + 2hiwo — 0. (30) 

Assuming w(^) 7^ 0, we normalize w{z) by wq = 1. Then, 
at given wi and W2 the remaining coefficients w„ are 
determined from Eqs. H29|l . The choice of wi and W2 is 
not arbitrary but determined by Eq. (|25|l . Thus, setting 
C = Ci and, hence, go — 0, we must choose wi in such 
a way that the series for w (or its continuation) 
converges at the second turning point C21 or satisfies the 
boundary condition at infinity l|25j) if C2 = +00. For the 
Mathieu and Lame equations this leads to the function 
w{z) which is an integral one, for the latter it being a 
polinomial 36j . In these cases the domain D is evidently 
the z plane with \z\ < 00. 

Let us suppose the function wiz) (|27|l is found. Return 
now to Eqs. (pni) - (|^ . Denote as W the Wronskian of 
the solutions (|^ . 



X+X^t - X+tX^ ^ const . 



(31) 



Setting 

yi = X, 
yi = x^ 

we then obtain 



7/2 = (zt^O), 

y2 = X+ [zt 0), (32) 



2/12/2 - 2/i2/2 = W/y^, 



(33) 



wher e ^ ^ is assumed. The system of equations H27f) 
and H33|l can be easily solved, which gives 



2/1,2 



exp 



wy/g 



: dz. 



where 



Now let us insert yi^2 back into Eq. it^ . 



(34) 

(35) 
We 



obta 



2gww" + g'ww' - gw'^ + Ahw^ + ^Q. (36) 

By this formula one can find the constant W'^ from a 
knowledge of ?x;(2;) in a vicinity of any point z. Thus, cal- 
culating H36|) at any one finite turning point C we obtain 
(in terms of expansions (|28|l with normalization i^o = 1) 

= -Aho - giwi. (37) 

Alternatively, one can take zeros Zi of u'(z). (The func- 
tions g{z) and 'w{z) do not have common zeros because 
otherwise the Wronskian (|31l) would be zero.) Then we 
find 

=g{zi)w'^{z,). (38) 

The requirement for positivity of determines the val- 
ues of parameters of Eq. (|20(l (resonance zones) for which 
the resonant solutions exist. 

Let us construct these solutions. Consider the intervals 
of monotonicity ti ^ t ^ t2 [zt ^ 0), t2 ^ t ^ ts {zt ^ 0), 
etc. According to (|32ll and l|34() we can write 

h 



Xi 



^±(ii)cxp 



Ci 



Xlit) - Xl{t2)exp 



C2 



<^V9 
f± 



dz, 



dz, etc. 



(39) 



To find the characteristic exponent n consider, e.g., the 
growing solution X_|_(/). Setting i = ^2, -^(^2) = C2 the 
first equation of and t = t^, 2; (is) = Ci in the second 

Using E 
-r, (p{t + T) 



one, we can express ^^(ts) through X'^(ti). Using Eq. 



(I21|l and taking into account that is 
±ip(t), we thus obtain 



ti+' 



dz 



(40) 



Recall that r is the period of z{t), the constant W is 
determined from Eq. (p?7jl or Eq. its sign being 

taken opposite to that of the integral in H40|l to provide 
for positivity of /i. Since w{z) has zeros, the integrals 
in Eqs. H39|l and H40(l are understood as their principal 
values. Formula H40|l is a sirnple generalization of the 
ones used previously in Refs. Issl [39i. t40j . 
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B. Evaluation of the characteristic exponent 

Let us return to Eq. H18|). It can be written in the 
form of Eq. 12(J|) if we set 



zit) 
h{z) 



- ln(a/a,] 



z(0) = 0, 



- 3 + + z. 

Equation Q then immediately gives 

g{z) = A[ul{e^ - 1) - z] 



(41) 
(42) 



(43) 



Zeros of this function are shown in Fig. 3. Since 
^(i) = a/flniax oscillates with the period T in the in- 
terval -1 < ^ < 1 [see Eqs. Q-ini)], the function 
z{t) (|41l) oscillates with the period t = T/2 between the 
turning points = and C2 ^ +00. To calculate the 

Im z 



Rez 



FIG. 3: The layout of zeros of g{z) 1431 1 on the complex z 
plane. The vicinity D (shaded) of the interval [0, 00) belongs 
to the domain of analyticity of wlz). 

characteristic exponent by the formula (|40|l we need to 
know the function w{z) which is the solution of Eq. H2t)|) 
with boundary conditions (|25|l . Unfortunately, for given 
h{z) (|42|l and g{z) (|43() equation (|26|l cannot be solved 
analytically. We solve it numerically for various values 



of the parameters E and = 1 — In a: 



2 

max 



Doing 



so, we use the conditions (|25|l in the following way. As 
discussed above, the fulfilment of 125() at a finite turning 
point means analyticity of w{z) in some vicinity of this 
point. Therefore, we can use the expansions setting 
there C = Ci = 0, .go = 0, gi = A{ujI - 1), = Aujl/n\ 



(n = 2,3,. 
then gives 



.),hQ^E + i 



3, hi = 1. Equation ^ 



W2 



1 + {2E + 3w§ - 6)wi 
6(^f5 - 1) ■ 



(44) 



We thus solve Eq. H26() with the following conditions at 
z = 0: w(0) = 1, w'{Q) ~ wi, w"(0) = 2^2. Given values 
of E and ujq, we choose wi so as to satisfy the condition 
at infinity, 



(45) 



At the same time, since fj, assumed to be real, the values 
of E, Wg, and Wi must provide for positivity of W^, 

^ A[3 - E - ujI - {luI - l)wi] > (46) 

[see Eqs. ((33 and (gOJ)]. Conditions (gSJl and igHI) deter- 
mine the resonance zones in the space of parameters E 
and luq (or a^^^). Hereinafter the zones will be referred 
to as Zj and numbered sequentially as E grows (with 
flj^jj^ fixed) starting with j = — 1 in the region E < 0. 
Figure 4 shows the solutions w{z) for zones Zi, Z2, and 
Z3 lying in the region E > 2. Now, knowing w{z), we can 




FIG. 4: 



onance zones 
2 

'^max 



Behavior of w{z) for E 
(a) zone Zi, 



2 

^max 

0.9025, (c) zone Z3, af-^^-^ 



4 in the different res- 
= 0.5625, (b) zone Z2, 
= 0.9598. The values 
of f^max choosen correspond to the centers of zones where fj, 
achieves its maxima. One can see that the number of zeros of 
w{z) is unit above the number of a zone. 

calculate the integral in (QDJ. Because w{z) has zeros, we 
first transform the integrand with the help of Eq. 136|) 
extracting the total derivative. Owing to the condition 
(|45|l , the latter does not contribute to the principal value 
of the integral, while the remaining terms give 



dz 



2W^ 



(w' y^)' \nw'^ + 8hw 



dz 



V9 
(47) 

The integrand on the right-hand side of Eq. (|47|l is more 
convinient for numerical integration because its singular- 
ities are all integrable. We perform the integration in 
(|T7jl . calculate by the formula and find the pe- 
riod T = 2r by integration of Eq. {Tj) . These procedures 
are carried out numerically for a set of grid points in ev- 
ery resonance zone. In this way from H4U|) we obtain the 
characteristic exponent /i as a function of E and a^^^. 

To check this result we derive ii{E, a^j^^) directly from 
analysis of numerical solutions of Eq. (|18|l . Examples of 
these solutions for resonance zones Zi, Z2, and Z3 are 
shown in Fig. 5. The growth of the amplitude with time 
is clearly seen. The function ii{E, a^^^^) so derived is 
found to be fully coincident with the one obtained by the 
formula 




t 



FIG. 5: Resonant solutions of Eq. (EHJ, X{t) ~ X+{t) = 
ip{t)e'^*, (solid lines) and the function z{t) — — ln(a/amax)^ 
(dashed lines): (a) zone Zi, (b) zone Z2, (c) zone Z3. The 
initial conditions are: a(0) = flmax, at(0) = 0, Xt{0) — 0, 
X(t) is normalized in a proper way. The values of E and 
flmax in each zone are the same as in Fig 4. It is seen that 
ifi{t) is r-periodic in Zi, r-antiperiodic in Z2, r-periodic in 
Z3, and so on, in accordance with the solutions 13911 [see Eqs. 
and ^ and Fig. 4]. 



The resulting stability-instability chart is presented in 
Fig. 6. Figure 6(a) shows the region E > 2. There is 
an infinite series of narrow resonance zones Zi, Z2, Z3, 
. . . , the first one having the highest magnitude (« 0.08 
at the maximum). All these zones originate from the 
point E ^ 2, al^^^ = 1 at which ^ = [see Eqs. (gO)) 
and (|46|l ] . In the region E ^ 2 we have two zones, Zq 
and Z^i, lying in the ranges < i? < 2 and E < 0, 
correspondingly. Since in these zones the values of fi 
proved to be much greater than in , .Z2 , . . . , we depict 



FIG. 6: The stability-instability chart: 
zones are shown, (b) E ^ 2. 



(a) E > 2, first ten 



C. Spatial structure of the perturbation 

Consider now Eq. (|19l) . It has the form of the 
Schrodinger equation for a quantum particle of the en- 
ergy E moving in the potential — Inu^. Since the poten- 
tial tends to -I-cxd with growing r [as r^, see Eqs. (|I3f) and 
(fnji ]. the energy spectrum is discrete, E = En, and the 
corresponding eigenfunctions ^n(r) are all localized. In 
the case of the nodeless pulson we have the isotropic 
harmonic oscillator. Its eigenfunctions are well known 
(see, e.g., 01 )• We write them as follows: 



vl'„(r) = Y.^l + i-ir-']R,u{r)Yii0,ip), (48) 
Rni{r) 



r^e'^'/Hi-^, l + l,r'], (49) 



P;™ (cos6')f 



(50) 



Here $(0,7,0:) is the Kummer function, PH'ix) are 
the associated Legendre functions, ci^m are constants, 
Ci,-m = (^i m- energy spectrum is given by 



E = E„ ^ 2n 



(n-0,1,2,...). 



(51) 



the surface ^J-{E, a^^^) for this region separately, in Fig. (Our energy levels are shifted with respect to the conven- 



6(b). 



tional ones since the minimum of the potential — In Uq is 
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-3.) 

In the case of the nodal pulsons the picture becomes 
more compUcated due to the loss of the orbital degener- 
acy. The corresponding eigenfunctions and eigenvalues 
can be calculated only numerically. As an example, in 
Fig. 7 is shown the energy spectrum for perturbations of 
the one-nodal pulson. 
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FIG. 7: The energy levels for the perturbations with different 
orbital numbers I. The case of the one-nodal pulson. 

Note that there always exist the eigenvalues E — 0, 
I = and the corresponding eigenfunction '9o(j') oc u{r). 
This fact immediately follows from the comparison of 
Eqs. lO and l(T^ . The corresponding Xo{t) in rj(t,v) 
is an oscillating function with the amplitude growing lin- 
early with time. It is easy to see, however, that this mode 
is physically meaningless. Indeed, it will formally appear 
if we perturb the pulson by a small variation of its am- 
plitude Omax but not the form u{r). Due to nonlinearity, 
this results in a pulson with slightly shifted frequency. 
Then the difference of the perturbed and unperturbed 
pulsons, i.e., ?7(i, r), will have the form of beats generated 
by two oscillations with close frequencies and the same 
profile u{r). The function XQ{t) approximates the initial, 
linearly growing part of a beat. We exclude this mode 
from the subsequent consideration, since it belongs to the 
class of perturbations that conserve a pulson as a whole. 
Next, for the nodal pulsons only, there is a mode with 
E = I = 1 (see Fig. 7). Since this mode cannot grow 
faster than linearly in time, we also do not take it into 
account. Further, we should exclude the mode resulting 
from a small translation of the pulson. The correspond- 
ing eigenfunction is proportional to nVu, where n is a 
displacement vector. Using Eqs. © and H19|) one can 
easily show that this mode corresponds to _E = 2, / = 1. 
Thus the resulting perturbation is written as 

r,(i,r) =^X„(t)\I'„(r), (52) 

n 

where X„ is a solution of Eq. (|18|l with E — £'„, _E„ ^ 
0, 2. If En and a^^^ are in a resonance zone, Xn{t) will 
be represented as a linear combination of the solutions 
(|21ll and, hence, will grow with time as e***^^" 



D. Instability of the pulsons 

The arrangement of the resonance zones on the 
(-B, a^j^x) plane indicates that for any spectrum i?„ there 
always exist the ranges of a^^x where pulsons are unsta- 
ble. But do the values of a,^„ax exist for which the pulsons 
are stable? To answer this question let us return to the 
surface ^Ji{E, a^^^ depicted in Fig. 6. Take, at first, the 
spectrum for the nodeless pulson. We choose the sec- 
tions /i„(amax) of the surface ijl(E, a^^^) hy E = 2n and 
project them on the (/.t, a^ax) plane. As a result, the 
pattern shown in Fig. 8(a) emerges. It is clearly seen the 




0.2 0.4 0.6 0.8 1.0 



2 



/,< 1 








0.2 0.4 0.6 0.8 1.0 



max 

FIG. 8; Superposition of the sections /i„(fflmax) of the surface 
fJ-{E,a'^i^^): (a) the nodeless pulson, E„ =4, 6, 8, 50 
(n = 2, 3, 4, . . . , 25), (b) the one-nodal pulson. En = -0.7142, 
0.4833, 1.1222, 1.8996, .... 

tendency to the total filling of the interval < a^^^ < 1 
by the resonant peaks as the successively higher en- 
ergy levels are accounted for. This implies that for any 
given a^^^ there always exists an unstable mode with 
/i = /i„(amax): i'^-, strictly speaking, all nodeless pul- 
sons are unstable. On the other hand, the figure shows 
that there are domains of a^^^ where the peaks are very 
small. These domains are the gaps between the main 
peaks originated from the low-energy cross sections of 
the surface ^{E, a^ax) over a few first zones. In the gaps 
the exponent /i is small, so that the corresponding pul- 
sons are long-lived. For example, in Ref. [33 we observed 
numerically the nodeless pulson with a^^x = 0.49 that 
conserved its coherency against the radially symmetric 
perturbations over the course of several hundreds of pe- 
riods. 

Further, the above projective procedure is performed 
using the spectrum of the one- nodal pulson (Fig. 7) . The 
main contribution here is made by the sections with the 
energies En = -0.7142 {I = 0), 0.4833 (/ = 2), 1.1222 
{I = 3), and 1.8996 {I = 4) falling into zones Z_i and Zq, 
the projections of the first and the third sections over- 
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lapping the other ones. The result is presented in Fig. 
8(b). We see that totally full. Thus, the 

one-nodal pulson has neither stability nor even quasista- 
bility domains. It seems likely that things will get worse, 
not better, if one goes to the multinodal pulsons. We 
thus conclude that, strictly speaking, all pulsons of the 
model considered are unstable. But nodeless pulsons can 
be quasistable in narrow ranges of amplitudes. It is the 
long-lived pulsons that can be of astrophysical and cos- 
mological interest. If the dark matter consists of scalar 
particles, such pulsons will be realistic candidates for the 
dark matter objects having oscillating density |4^ . 

E. On the instability of the Affleck-Dine type 
condensate 

The obtained stability-instability chart turns out to 
be appropriate for the stability analysis of the nonlocal- 
ized coherent states as well. As an example, we consider 
a uniformly distributed background 0o(Oj * scalar con- 
densate, oscillating around the minimum of the poten- 
tial Q at (/) = 0. This state can be formally obtained 
from Eq. |^ if we set there u{r) = 1. We thus assume 
that </>o(i) obeys Eq. Taking the perturbed state 

4> = 00 (^) I") 5 ill the linear approximation from Eqs. 
(0) and l@J we readily obtain 

Att + ik^ -2-lncj)l)A^0, (53) 

where A{t,'k) is the Fourier amplitude of the perturba- 
tion, and k = |k|. It is seen that the real and imagi- 
nary parts of this equation have the form of Eq. H2()|l 
with h{z) given by Eq. z = , 

cJq = 1 — In^Qjjjj^^, and E ~ k'^ . Returning to the 
stability-instability chart (Fig. 6) we note that in the 
region E ^ maximal values of /i are attained in the 
zone Zq for which < i? < 2. Interestingly, this band 
exactly coincides with the one obtained in Ref. for 
the power-law potential approximating when A 1. 
In the interior of Zq the exponent n depends almost not 
at all on the amplitude of the condensate oscillations and 
is a sufficiently smooth function of fc^ with a maximum 
at fc^ = fco ~ 1 where fi w 0.5. Therefore, if the initial 
power spectrum \A{0, k)| lies in the region < |k| < \/2 
and, in addition, its characteristic width along k is small, 
Afc <^ v^j then the growth of the perturbation ampli- 
tude will not be accompanied by significant changes in 
the structure of the perturbation. The limiting case of 
such perturbations is a harmonic wave. Otherwise, if 
Afc > \/2, the shape of the power spectrum will vary 
with time so that a maximum will appear at fco ~ 1. As 
a result, the effective width of the spectrum will become 
smaller, Afc < 1. In this case, if the initial spectrum is 
sufficiently isotropic in k space, the parametric amplifi- 
cation of the perturbations will result in the emergence of 
the localized field configurations of the characteristic size 
Ar ^ 1/ Afc > 1 that agrees with the radius of the gaus- 
son (see Sec. II). At this scale the field practically does 



not undergo spatial oscillations since the corresponding 
wavelength 27r/fco > 1. We thus expect that at the non- 
linear stage these configurations will turn into the node- 
less pulsons. Their period will be equal to the period 
of the condensate oscillations since in the zone Zq the 
parametric amplification proceeds at the basic frequency. 
Gradually, the energy of the oscillating condensate will 
go to ensemble of the arising pulsons, this process result- 
ing in the damping of the background oscillations. As to 
the pulsons themselves, they can be long-lived or short- 
lived depending on their amplitudes, in accordance with 
the results of the previous Subsection. 

Note, that numerical simulations performed for the 
complex version of the model (0) and JHl have shown 
the fragmentation of both the rotating 18, 1^ and os- 
cillating i2Qj Affieck-Dine condensate. The localized con- 
figurations arising in the condensate have been identified 
with Q-balls. We belive, however, that the configura- 
tions observed in the oscillating condensate are in fact 
the complex pulsons (see Sec. IV), rather than the usual 
Q-balls. This possibility was early discussed in Ref. 
where an attempt to simulate the complex pulson has 
been made. 

The resonant exitation of the pulsons was also observed 
in the two-vacuum (j)'^ — (p^ model within a regularly os- 
cillating background 39] and in the model within an 
initially thcrmalized background ji^. Note that in two- 
vacuum models the pulsons can play the role of nuclei 
of a new phase. In Ref. the general suggestion has 
been made that the parametric resonance can underlie 
the mechanism responsible for the first-order phase tran- 
sitions in nonlinear non-dissipative systems. This con- 
jecture turns out to be in agreement with recent results 
of Ref. ^3 where the resonant nucleation within the 
thermalized background have been numerically observed 
in the (f)'^ — 0^ model. Note, in addition, that the dy- 
namical nucleation can also take place in the nonlinear 
Schrodinger equation [45l| . 



IV. CONCLUDING REMARKS 

In this paper we have examined only the linear stage 
of instability at which small deformations of the pulson's 
shape result in loss of the coherence. There is numerical 
evidence that in time the growth of the perturbations 
becomes saturated due to nonlinear effects . We thus 
suggest that in the model considered the pulsons, while 
unstable, remain well localized objects with no tendency 
for spreading or collapsing. 

Further, we dealt with a real scalar field. It would be 
interesting to perform the similar analysis for a complex 
scalar field too. It is believed that the existence of the 
scalar charge can stabilize a field lump. For CKballs this 
fact is well established (so-called Q-theorem QjBlls)- 
contrast, for the complex pulsons this is an open question. 
As it was shown in Refs. ^,^2§\, the field equation Q 
with U' — — 01n(00*) admits the exact pulson solutions 
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of the form cj)Q{t,r) = a{t)u{r)e^^^*\ where a{t), u{r), 
and 9{t) are reaL The function u{r) satisfies Eq. © as 
before, while a{t) oscillates with a period T in accordance 
with the equation 



att = -3- 
aa 



(54) 



where 5 is a real constant, < (2e)~^, and 9t = qa^^. 
The constant q is proportional to the charge of the scalar 
field. In contrast to Eq. ||SJ|, the potential in the square 
brackets of Eq. H54I) prevents a{t) from being zero. With- 
out loss of generality one may assume a{t) positive, so 
that the oscillations occur arround the minimum of the 
potential at a = ao, where oq is the least positive root 
of the equation a^lna'^ — ~q^. If a is at rest in this 
minimum, then 9{t) = qa^^t + 9(0), and we have the 
standard Q-ball. Physically, Eq. H54I) describes the mo- 
tion of a mechanical particle with an angular momentum 
q in the potential (a^/2)(l-lna^) The condition for 
its trajectory to be closed is 6{T) — 9{0) = 2TTm/n, where 
m and n are arbitrary integers. In fact, it relates the en- 
ergy of the particle and its angular momentum whereby 
such trajectories exist. In our case this means periodicity 
of the solution (1)^(1, r) with the period nT. Obviously, 
there is an infinity of such solutions. Taking (j)Q{t,r) and 
considering the partial perturbation 77 a X(t)^(r) one 
can find that the function Vl/(r), assumed to be real, sat- 
isfies Eq. H19|l as before, while X{t) obeys the equation 



{E-l~liia^)X = e''"'X 



(55) 



where E is a real constant. This equation can be rep- 
resented as a system of four real first-order equations 
with periodic coefficients of the periods T and nT. It 



is significant that, since a{t) ^ 0, these coefficients are 
bounded in time, so that one can attempt to estimate the 
characteristic exponent of the system using the standard 
methods |47j . 

Also, it would be interesting to examine stability of a 
selfgravitating pulson. Hopefully, gravitation can expand 
the domains of (quasi)stability, as it is the case for Q-balls 
4]. These are possible subjects of our future work. 

In the present paper we have investigated stability of 
both the coherent localized states (pulsons) and nonlo- 
calized states (uniformly oscillating scalar condensate) 
of the real scalar field. Our main analytical result is the 
generalization of the Lindemann-Stieltjes method to the 
case that the periodic coefficient in the Hill equation is 
unbounded in time. Our main numerical result is the 
stability-instability chart with the values of characteris- 
tic exponent calculated in the resonance zones. Using 
this chart we have found the gaps in the set of the pul- 
son amplitude values in which the real nodeless pulsons 
conserve the coherency for an extremely long time. Also, 
considering the oscillating scalar condensate, we have de- 
temined the wavelength of the most unstable mode. This 
wavelength turned out to be equal to the characteristic 
size of the nodeless pulson. We thus suggest the pul- 
sons can be formed due to resonant fragmentation of the 
scalar condensate. These are our main physical results. 
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